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Abstract 


In Argentinean Patagonia, there exists a growing interest in understanding local problems 
associated with the increasing of wind energy penetration in the grid. Recent papers show that some 
of these power quality problems can be reduced by the proper control of the wind farms. In this way, 
this work deals with the impact and potential contribution of wind generation on the damping of the 
electromechanical oscillations called inter- and intra-area oscillations. For gaining qualitative 
insights and understandings on this complex subject, a test system that allows to isolate the 
oscillations modes is considered. The analysis shows that new control concepts for wind farm can 
efficiently contribute to the power system damping. 
© 2007 Elsevier Ltd. All rights reserved. 
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1. Introduction 


From early times many countries have been modifying their legislation in order to 
promote the development of renewable energy. Then, political and economical incentives 
tend to encourage active power generation as much as possible from wind farms. As 
examples, green certificates, different kinds of tax exemptions and financial helps to the 
wind power generation are applied in Chubut State (Argentine) where the biggest 
argentinean wind farm is located. 

However, as the number and size of wind farms increase, the interaction between wind 
farms and grids becomes more important [1—4]. Then, associated with the total amount of 
power connected to a distribution network, some power quality problems can arise due to 
the fluctuating nature of the wind, e.g. voltage and power disturbances which deteriorates 
the network performance. This fact, which has limited the wind farm installations, is 
potentiated in weak grids in the southern tip of South America in Patagonia, where the 
electrical network is not connected to the main interconnected argentinean power system 
and where the medium wind velocities can be over 12 m/s [5]. 

These characteristics have generated in Patagonia the necessity of studying the impact of 
wind generation on electrical networks and researching about the potential contribution of 
wind farms to overcome some network problems. In this way, this work resumes studies 
which show how the wind farm control could contribute in one of the aspects of power 
quality: the damping of electromechanical oscillations. 

Oscillations are undesirable because they limit power transfers on transmission lines 
and, in some cases, induce stress in the mechanical shaft of synchronous generators [6]. 
These oscillations are of two types [6,7]: 


è Local mode or intra-area oscillations. They are associated with rotor angle oscillations, 
of a single generator or a single plant against the rest of the power system. They usually 
have frequencies from 1 to 3 Hz. 
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èe Global mode or inter-area oscillations. They involve oscillations of a group of 
generators in one area against a group of generators in another area. Inter-area 
oscillations are in the range of less than 1 Hz. 


In conventional networks (i.e. without renewable energy resources contribution), small 
signal analysis and control are well-established techniques [6-8]. These techniques are 
employed in this work looking for obtaining indexes which allow to evaluate the potential 
contribution of wind farms for damping electromechanical oscillations. The analysis 
considers different percentages of wind power penetrations and different kinds of basic 
control strategies. Simulations are presented to support the analytical results. 

The structure of this paper is as follows. In Section 2, a wind farm model is presented. In 
Section 3, the wind farm model is included in the reference frame of a power system model. 
Section 4 presents the wind farm control which considers different strategies, particularly it 
includes simple control strategies to regulate frequency and voltage. In Section 5, the 
impact of wind farms in a power system is evaluated through an eigenvalue study, and 
simulation results are carried out. Finally, in Section 6 conclusions are presented. 


2. Model of a wind farm 


A complete model of a wind farm with a high number of windmills may lead to compute 
an excessive and impractical number of equations. The size of the wind farm model may be 
reduced by aggregating several wind turbines with similar incoming wind into a bigger 
turbine called aggregated turbine [9]. The mechanical and electric parameters per unit are 
preserved, and the nominal power is increased up to the sum of the nominal power of the 
whole set of turbines to obtain the parameters of the aggregated turbine. 

Making a decision about how many turbines should be aggregated depends on each 
particular case. Analyzing network changes far from the wind farm, it could be enough to study 
the wind farm as an only one aggregated turbine. However, when a disturbance takes place in 
the wind farm it could be necessary the study of all the generators in a separate way [10]. 

As the main objective of this paper is to examine the effects that the wind farm has over 
the network oscillations, a single lumped mass of the aggregated wind turbine is considered 
[11-13]: 


Ji : Ji\; Tt 
(n+ i)n = (nti) T (1) 


which is referred to generator shaft by the gearbox transmission N, being J; and J, the 
inertias of the aggregated wind turbine and the generator plus the gearbox one, and Q, and 
Q, generator and turbine speeds, respectively. 7; is the turbine torque (Appendix A): 


npr? 
T2 
and Ts is the generator torque defined in Appendix B. Because of in Patagonia most of the 
wind generators are powered by double fed induction (wound rotor) generators (DFIGs), 
in this work these electric machines are the only ones considered, even when the 
conclusions can be applied to other variable speed generators. Due to the electric dynamics 
is fast against that one of the rest of the system [14,15], the dynamical aggregated model of 
a wind farm can be approximated by the corresponding mechanical behavior. 


Tı wCp(A)/Qt, (2) 
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3. Power system model 


The linearized differentialalgebraic (DAE) model for a conventional multimachine 
power system (without wind farms) is [6] 


Ax = A\Ax + BAI, + BAV; + E\Au, (3) 
0 = C)Ax + Di Al, + Dr AV zg, (4) 
0 = CpAx + D3AIg + D4AV g + DAV, (5) 
0 = DAV; + DAN, (6) 


where the expressions (3)-(6) represent synchronous generator dynamics, static stator 
synchronous generator equations, active and reactive power equations at generator buses 
and at load ones, respectively; AT = Lan deal are generator currents; AV, = [0n Vale 
are phases and amplitudes of voltages at generator buses; Au=[w; Tmi Vrn]! are the 
synchronous frequency, mechanical torques and reference voltages, and AV; = [0n Vn]' 
are phases and amplitudes of voltages at nongenerator buses. 

In order to include the wind farm in the electrical network equations (3)-(6), the 
linearized model of the aggregated DFIG wind farm operated at constant wind speed is 
calculated from expression (1): 


JAi\.. AT 2H 
@ + z) AQ, = 1 L AT, => E Ao, 
N Mrnom 
— (Ktur/N?)(2/P)Ao, AT: 


7 PBASE /Ornom(2/P) PBASE /Ormnom(2/P) l u) 


Eq. (7) is expressed in p.u., being Hyg = 1/2 (Jg + (J /N7))(@rnom(2/P))”)/ Pass, co, the 
electrical generator speed, W;nom the rated electrical generator speed and Pgase the rated 
power for scaling all equations in per unit. Also, the gain Kurp is obtained from Appendix 
A, taking the derivative of Eq. (16): 


_ (npr? /2)v3AC,/AQ, — (apr? /2)? Cpi 


K turb OQ 
tl 


(8) 


where the subscript “‘l” indicates the linearization point. 

In order to model the complete system (with wind farms), it is necessary to modify 
expressions (5) (or (6) if the wind farm is connected at a nongenerator bus) with the 
generated active power and to modify expression (3) as a consequence of the new state, the 
wind farm equivalent speed (œr). Also, the generator torque takes different expressions 
depending on the control strategy applied. Then, before proceeding with the model (AT.), 
it is important to analyze the possibilities about the control of the active and reactive 
powers of the wind farm. 
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4. DFIG wind farm control strategies 


As it is pointed in Appendix B, vector control can drive independently active and 
reactive powers of DFIG. Then, this work considers the next cases: 


e CP: Active power control (tracking or regulation) and reactive power control AQ, = 0. 

e FC: Active power control APxas) (in order to contribute to frequency regulation) and 
reactive power control AQ, = 0. 

e FCVC: Active power control AP«ar) (in order to contribute to frequency regulation) 
and reactive power control AQ) (voltage control). 


4.1. Conventional wind farm strategies (CP) 


One of the most conventional modes of operation of wind farms powered with DFIGs 
has been to track the curve Tio (expression (20), Appendix A): 


Tl r? r 3 
j (=) Ca (9) 


thus getting the maximum active power available in steady state: Pio = Ti Qt. 

In this operation mode, wind farms disturb the network as the wind velocity changes. 
Then, given the growing contribution of wind generation in electrical systems the 
aforementioned strategy of maximizing the power conversion can be unsuitable for some 
cases. Hence, other strategies are considered such as the wind farm power regulation at a 
power reference (Psc) [16]. 

Obviously, no contribution to reject frequency disturbances can be expected from any 
two previous wind farm control strategies. This is because of the fast power (vector) 
control which keeps the wind farm power at its reference value. 


4.2. DFIG wind farms with frequency control (FC) 


In order to develop a corrective action to frequency changes, the next simple 
modification of the power reference in the DFIG wind farm is evaluated [17]: 


Pyt = Psc + APaas), (10) 
APaan = KAF, (11) 


being Pyr the wind farm power reference, Psc the static power reference, and AP«af) a 
variable component of the wind power reference which is proportional to the frequency 
variation Af. 

In order to assure the proposed control strategy, it is necessary to consider an inertia 
reserve of the wind farm from operating each wind turbine with a power coefficient Cp 
lesser than the optimum one (Appendix A). This fact implies that a static power reference 
Psc < Pio must be assigned at the nominal frequency f. Consequently, the wind generator 
has 100 * (Cyo — Cp1)/Cpo% of additional power to contribute to the frequency 
regulation in steady state. However, it its important to note that during a transient, the 
wind farm can contribute with the excess of energy available in the rotor inertia. 
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Fig. 1. Power—speed—frequency curves. 


Fig. 1 shows power-turbine speed—frequency curves with the wind velocity as a 
parameter. In particular, wind turbine characteristic curves (expression (17) in Appendix A) 
are depicted in the power-speed plane, and frequency characteristics are presented in the 
power-frequency plane. Note that power—frequency characteristics are equivalent to the 
speed droop characteristics in synchronous generators [8,7], however, the slope of 
power—frequency characteristics of the FC wind farm can be more pronounced. In this 
figure the operation point A corresponds to the Psc for nominal frequency f, and 
Cpi < Cpo. Because of the turbine speed corresponding to the point A’ is larger than the 
turbine speed in point B’, there exists an excess of energy stored into the turbine inertia. On 
the other hand, on each power—frequency characteristic, exist two upper limits, the static 
limit B” and the dynamic limit C”. The dynamic limit C”, in correspondence with Pio, is 
specified for preserving the wind farm operation when the frequency remains below its 
nominal value. Then, it is possible to use the energy in between two curves, the first curve is 
the turbine speed characteristic P, and the second one is due to Pig which correspond to 
Eqs. (17) and (19) in Appendix A, respectively. 

In the present work, the lower limit considered for the generated power is zero. Then, if 
frequency increases, the wind farm can reduce its active power injection into the network 
allowing that the energy can be stored in the inertia, then, all the excess of energy goes to 
the inertia reserve. Hence: 


0< Pwr < Pioo. (12) 


Contrary to the conventional synchronous generation, double fed induction machines 
are able to modify very quickly the generated power in comparison with the rest 
of the electrical system [15]. The aforementioned implies that the participation of 
the wind farm could be effective, even when the wind power involved would be 
significantly lesser (for reasons of power quality [10]) than the conventional installed 
power. 
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When a DFIG wind farm is operated with FC, expression (7) is modified with: 
Ate Ais (Kturp/N?)(2/P)A@, Kto, AO: + Kto, A0; 
@rnom i PBASE /Ornom(2/P) PBASE /Ornom(2/P) Í 


being Kto, = —Psc/w2, and Kto, the FC gain. 
Then, for power analysis, expression (5) (or (6)) is modified adding the linearized expres- 
sion for the FC wind farm: 


APeafy = K Po, AOe, (14) 


(13) 


considering K po, = Kto, Qro. 
4.3. DFIG wind farms with frequency and voltage control (FCVC) 


About the control of the wind farm reactive power AQ,, strategies accepted nowadays 
are: 


è Q, regulation, i.e. AQ, = 0 (assumed in previous cases); 
è cos ọ regulation, i.e. cos ọ = specified [18]. 


However, taking into consideration the increasing penetration of wind farm generation 
and looking for satisfying the following objectives [7,19]: 


o keeping voltages at the terminals of the equipment within acceptable limits, 

o enhancing system stability to maximize utilization of the transmission system, 

o minimizing the reactive flow so as to reduce losses ensuring the transmission system 
operates mainly for active power transfer, 

e a DFIG wind farm voltage regulation is evaluated, i.e. AQ, = AQxy). 


Considering active compensation, devices as synchronous compensators and SVCs 
automatically adjust the reactive power absorbed to maintain voltages of the bus. From 
this viewpoint, DFIGs with reactive power control can also be considered as active 
compensation devices. 

Taking into account that the wind farm reactive power can be controlled independently 
of the active one, expression (5) (or (6)) is modified to take into account the voltage 
control. Here, the most simple case is considered: 


AQay) = KovAV, (15) 
being Kgy the control gain. 
Note that, as with synchronous compensators, the reactive power injection will be 
limited by the nominal apparent power S, of the DFIG [18]. 
5. Impact and potential contribution of a wind farm in a power system 


5.1. Selected test system 


In this section, it is evaluated, how basic control strategies of the DFIG can affect the 
dynamical behavior of a power system. In this way, Fig. 2 depicts the selected test system. 
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Fig. 2. System under test: inter- and intra-area power flows. Thick line: intra-area oscillation. Dashed line: inter- 
area oscillation. 


It is formed by two areas. Area 1 is an infinite bus meanwhile Area 2 is formed by two 
synchronous machines in parallel with a wind farm feeding a load which have been 
considered of constant impedance type. Each synchronous machine, where it is considered 
a constant motor torque, is associated with a standard IEEE Type I voltage regulator [7]. 

The symmetrical network topology considered in Fig. 2 is specially chosen for analyzing in 
a separate way intra- and inter-area oscillations [20]. This topology permits to simplify the 
analysis by separating the calculus of the network eigenvalues for each oscillation mode. 

In Fig. 2, it is considered that the disturbances appear in the loads at buses 1 and/or 1’. 
It is well known that small signal analysis allows splitting any disturbance input in 
different components and then to calculate the system response applying the superposition 
principle [21]. To the ends of the present study it is interesting to split the disturbance input in 
such a way that: 


@ one of this new components (named signal or disturbance or fault of common mode 
[22]) only excites the inter-area oscillation and 

e the other component (named signal or disturbance or fault of differential mode [22]) 
only excites the intra-area oscillation. 


In this way, and taking into account the network topology of Fig. 2, any fault at bus 1 and/ 
or 1’ (Py; and Py, respectively) will be divided in a common load signal [22]: 


Puy t+ Pry 


Pic= 7 3 
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which is applied at the same time at both buses 1 and 1’, and a differential load signal: [22] 


Pii — Puy 
Pip = L Li’ 


applied at the same time in both buses 1 and 1’ but with opposite signs in order to 
reconstruct Py; and Pry: 


Pur = Pic + Pip, 
Pir = Pic — Pip. 


Note that, in Fig. 2, inter-area oscillations, and no intra-area ones, are obtained from Pic 
(i.e. the component of fault which is equivalent to two equal signals with the same sign at 
buses 1 and 1’). In this case both generating groups (synchronous generators and wind 
farms) behave in the same way respect to the infinite bus, i.e. the synchronous generator 
shafts experience the same motions indicating that all the Area 2 behaves coherently 
against the Area 1. Then, the power flow in between the areas corresponds to the indicated 
in Fig. 2 with dashed lines. As a consequence of the network symmetry and the 
characteristic of the Pic component of the disturbance, there is no flow in between the 
synchronous generators (buses | and 1’), i.e. there are not intra-area oscillations. 

On the other hand, intra-area oscillations can be studied from the differential 
component of the disturbance which perturbs the network at buses 1 and 1’, the same 
quantity and at the same time but with opposite signs. Due to the characteristic of Pip and 
the network topology, the excess of power in one generator group in Area 2 is exactly 
absorbed by the lack of power on the other generator group in the same area. As a 
consequence of the power flow, which is depicted in thick line in Fig. 2, there is no 
variation of the power flow generated/absorbed by the infinite bus, i.e. there are not inter- 
area oscillations. 

Finally, the real response of the system is obtained, by applying the superposition 
principle, from the addition of the system responses to both components of the 
disturbance. 


5.2. Cases of study. Simulation results 


In order to evaluate the power system dynamic behavior when wind power penetration 
is increased (the same quantity that synchronous generation is decreased), three cases are 
studied by means of an eigenvalue analysis based on Fig. 2: 


1. wind farm with constant active and reactive power regulation (CP), i.e. AP. = 0 and 
AQ. = 0, 

2. wind farm with FC, i.e. AP. = AP«af) and AQ, = 0, 

3. wind farm with frequency and voltage control (FCVC), i.e. AP. = APyay) and 
AQ, = AQ. y)- 


Fig. 3 presents simulation results for the CP case taking into consideration the selected 
system in Fig. 2 with a 150MW load at both buses (1 and 1’) and wind farm and 
synchronous generations of 30 and 120 MW at each bus, respectively. A disturbance of 
40 MW at bus 1 is considered at t = 2s. The transient response of frequency at bus 1 is 
presented in part (a), this response is separated in intra- and inter-area oscillations in parts 
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Fig. 3. Dynamic response of the system considering the CP case and a wind farm penetration of 20%: (a) 
frequency at bus 1, (b) intra-area oscillation and (c) inter-area oscillation. 


(b) and (c), respectively. The intra-area oscillation corresponds to the difference between 
shaft velocities of synchronous generators at buses | and 1’ and the inter-area response is 
obtained by taking away the inter-area frequency of part (b) from the transient response of 
frequency of part (a). 

In next subsections, the effects of FC and FC and voltage control (FCVC) cases on both 
kinds of oscillations, intra- and inter-area ones are analyzed and, finally, the total response 
of the electromechanical oscillation, voltage and control effort are presented. 


5.3. Intra-area oscillations 


Based on the considerations commented in Section 5.1, the eigenvalues corresponding to 
the intra-area oscillations are calculated and presented in Table 1. Parts (a)-(c) of Table 1 
corresponds to CP, FC and FCVC cases, respectively. Note that, looking for emphasizing 
the problems under study, the eigenvalues, which are scaled by 27 giving frequency in Hz, 
are calculated for a broad span of wind penetrations. 

From analyzing Table 1 it is observed that: 


è At the same wind power penetrations, in parts (a)-(c), the most significative differences 
correspond to the eigenvalues in between lines. 
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Table 1 


Intra-area oscillation 
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Peol 0 MW Peol 30 MW Peol 60 MW Peol 90 MW Peol 120 MW Peol 150 MW 
(a) Case 1: DFIG wind farm with constant power and Q, = 0 

— 157.62 — 157.62 — 157.62 —157.62 —157.62 —157.63 
—431 —4.29 —4.27 —4.25 —4.25 —4.25 
—2.95 —2.95 —2.94 —2.94 —2.95 —2.95 
—0.42 —0.40 —0.39 —0.39 —0.39 —0.39 
+1.15i +1.15i +1.15i +1.14i +1.14i +1.14i 
—0.075 —0.093 —0.108 —0.117 —0.122 —0.123 
+1.42i +1.42i +1.42i +1.42i +1.43i +1.42i 

0 —0.18 —0.18 —0.18 —0.18 —0.18 
(b) Case 2: DFIG wind farm with frequency control and Q, = 0 

—157.62 —157.62 —157.62 —157.62 —157.62 — 157.63 
—431 —4.34 —4.51 —4.97 —5.96 —6.01 
—2.95 —2.95 —2.95 —2.94 —2.95 —2.95 
—0.42 —0.41 —0.44 —0.46 —0.43 —0.42 
+1.15i +1.14i +1.09i +1.2li +1.18i +1.17i 
—0.075 —0.17 —0.38 —0.70 —1.05 —1.05 
+1.42i +1.42i +1.38i +1.01i +1.49i +1.50i 
0 —0.18 —0.18 —0.18 —0.18 —0.18 
(c) Case 3: DFIG wind farm with frequency and voltage control 

— 157.62 — 157.62 — 157.62 —157.62 —157.62 — 157.63 
—4.31 —4.36 —4.55 —5.02 —6.01 —6.02 
—2.95 —3.12 —3.28 —3.35 —3.38 —3.38 
—0.42 —0.40 —0.43 —0.45 —0.39 —0.39 
+1.15i +1.13i +1.06i +1.2li +1.05i +1.04i 
—0.075 —0.17 —0.39 —0.70 —1.05 —1.06 
+1.42i +1.42i +1.371 +1.03i +1.49i +1.50i 
0 —0.18 —0.18 —0.18 —0.18 —0.18 


Eigenvalue analysis. 


e These eigenvalues present, at each wind power penetration, a higher damping factor for 
the FC and FCVC cases with respect to the CP one. This implies a potential 
improvement in the damping of the intra-area oscillation for the FC and FCVC cases. 

e For each control strategy, the damping of these eigenvalues increases with the wind 
penetration. This implies that the damping of the intra-area oscillation increases with wind 
penetration. Again, this increment is much more marked in the FC and FCVC cases. 

e In the CP case, the eigenvalues displacement with the increment of the wind energy is 
not significative. This fact is attributable to the independence of the CP strategy with 
respect to the frequency and voltage variations. In this case, the eigenvalues movement 
is a consequence of unloading the synchronous generators with wind penetration, i.e. 
there is a drop in the load which “sees” each synchronous generator provided that more 


wind energy is injected into the network. 
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Figs. 4 and 5 show the movement, for different wind energy penetrations, of the 
aforementioned eigenvalues for the CP and FC-FCVC cases, respectively. Fig. 5 indicates 
the strong (positive) influence on the damping of the intra-area oscillation with the FC and 
FCVC control strategies. 

Fig. 6 presents, for the three strategies, the intra-area oscillation when takes place a load 
change in 40 MW at t= 2s at bus 1, being the wind farm power generation and the 
synchronous power generation of 30 MW and 120MW at each bus, respectively. 
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Fig. 4. Intra-area eigenvalues for the CP case at different wind farm penetrations. 
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Fig. 5. Intra-area eigenvalues for FC and FCVC cases at different wind farm penetrations. 
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Fig. 6. Intra-area oscillation and voltage at bus 1 with a wind farm penetration of 20%. Wind farm with: constant 
power (CP), frequency control (FC) and frequency and voltage control (FCVC) strategies. 


Simulation results agree with the expected improvements from Table 1 and show the 
strong impact on the intra-area oscillation of the FC and FCVC control strategies when 
they are compared with the CP one. At the same time it is observed that the transient 
behavior of the intra-area mode is almost the same for the FC and FCVC control strategies 
but the last one, as a consequence of the reactive power loop (AQ.;)), produces a better 
voltage profile. 


5.4. Inter-area oscillations 


Table 2 presents, based on the considerations commented in Section 5.1, the eigenvalues 
corresponding to the inter-area oscillations for the three aforementioned cases. 
Analyzing Table 2 the next conclusions (similar to the Table 1) can be obtained: 


e Comparing the inter-area eigenvalues for the three control strategies at equal wind 
power penetrations, the main differences correspond to the eigenvalues in between lines. 

e Also, for each case, the damping of these eigenvalues increases with the wind 
penetration and it is much more marked in the FC and FCVC cases. This is shown in 
Figs. 7 and 8 where are presented the evolutions of the aforementioned eigenvalues for 
the CP and FC-FCVC cases, respectively. This implies that the damping of the inter- 
area oscillation increases with wind penetration. 

e These eigenvalues present, at each wind power penetration, a damping factor which is highly 
superior for the FC and FCVC cases with respect to the CP one. This implies a substantial 
improvement in the damping of the inter-area oscillation for the FC and FCVC cases. 
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Table 2 
Inter-area oscillation 


Peol 0 MW Peol 30 MW Peol 60 MW Peol 90 MW Peol 120 MW Peol 150 MW 


(a) Case 1: DFIG wind farm with constant power and Q, = 0 


—157.62 —157.62 —157.62 —157.62 —157.62 — 157.63 
—3.76 —3.71 —3.68 —3.66 —3.64 —3.64 
—2.70 —2.70 —2.70 —2.71 —2.72 —2.74 
—0.49 —0.48 —0.48 —0.48 —0.48 —0.49 
+1.46i +1.47i +1.48i +1.49i +1.50i +1.50i 
—0.039 —0.046 —0.051 —0.056 —0.059 —0.061 
+1.001 +1.00i +1.00i +0.99i +0.99i +0.99i 
0 —0.18 —0.18 —0.18 —0.18 —0.18 
(b) Case 2: DFIG wind farm with frequency control and Q, = 0 

—157.62 —157.62 —157.63 —157.63 —157.63 —157.63 
—3.76 —3.75 —3.90 —2.73 —2.38 —2.36 
—2.70 —2.70 —2.71 —1.60 —2.42 —2.41 
—0.49 —0.46 —0.46 —0.50 —0.51 —0.51 
+1.46i +1.48i +1.53i +1.53i +1.5li +1.50i 
—0.039 —0.20 —0.58 —4.49 —6.27 —6.38 
+1.00i +0.98i +0.74i &—1.60 &—0.23 &—0.23 
0 —0.18 —0.18 —0.18 —0.18 0.18 

(c) Case 3: DFIG wind farm with frequency and voltage control 

—157.62 —157.61 —157.61 —157.60 —157.60 —157.60 
—3.76 —3.77 —3.94 —3.34 —3.38 —3.38 
—2.70 —3.02 —3.26 —1.70 —2.53 —2.52 
—0.49 —0.46 —0.45 —0.47 —0.44 —0.44 
+1.46i +1.46i +1.47i +1.38i 1.271 +1.26i 
—0.039 —0.20 —0.58 —4.51 —6.26 —6.37 
+1.001 +0.98i +0.74i &—0.47 &—0.23 &—0.23 
0 —0.18 —0.18 —0.18 —0.18 —0.18 


Eigenvalue analysis. 


Simulation results, presented in Fig. 9, show the inter-area oscillation and the voltage at 
bus 1. A strong impact on the inter-area oscillation of the FC and FCVC control strategies 
when they are compared with the CP one is observed. Note that, the inter-area oscillation 
presents a slightly better behavior for the FC case when it is compared with the FCVC 
case. However, the last one improves the voltage profile. 


5.5. Complete transient response at bus 1 
Fig. 10 presents, for the three cases and at bus 1, in part (a) the complete transient 


response of the frequency, i.e. the shaft oscillation of the synchronous generator at bus 1, 
in part (b) the voltage and in part (c) the apparent power of the wind farm. 
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Fig. 7. Inter-area eigenvalues for the CP case at different wind farm penetrations. 
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Fig. 8. Inter-area eigenvalues for FC and FCVC cases at different wind farm penetrations. 


Parts (a) and (b) of Fig. 10, clearly show the improvements over the transient response of 
frequency and the voltage for the FC and FCVC cases when they are compared with the 
CP case. The reactive power loop (AQ,,;)) produces in part (a) of Fig. 10 a slight decrease 
in the performance of the transient response of the frequency with respect to the FC case. 
This fact is attributable to the lesser variation in the voltage observed in part (b). Indeed, as 
the voltage changes, constant impedance loads (considered in this work) tends to damp the 
frequency oscillation by modifying the active power absorbed. 

In part (c) of Fig. 10, it is noted that when the disturbance appears, DFIG wind farms 
with FC and FCVC react carrying out a contribution of active, and active and reactive 
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Fig. 9. Inter-area oscillation and voltage at bus 1 with a wind farm penetration of 20%. Wind farm with: constant 
power (CP), frequency control (FC) and frequency and voltage control (FCVC) strategies. 


powers, respectively. The reactive power contribution is corroborated in part (b) where is 
observed the improvement in the voltage profile. Note, that the maximum (transient) 
increment of the apparent power, i.e. the control effort involved, is about a fourth part of 
the steady state condition with less than a half second of duration. However, this is enough 
to strongly improve the transient response of frequency in part (a) and the voltage response 
of part (b). 


6. Conclusions 


In Argentinean Patagonia, there exists an increasing interest in understanding 
some problems associated with wind power generation. This fact is motivated for 
particular characteristics of the region as weak grids and high wind speeds. In this way, this 
work presents studies which show how control strategies of the wind farms could 
contribute in one of the aspects of power quality: the damping of electromechanical 
oscillations. Also, considering the apparent power available, a reactive power control 
strategy which enables to improve the voltage profile at the bus where the wind farm is 
connected was evaluated. 

Three cases were considered, wind farms with conventional power control (constant 
power: CP), with frequency control (FC) and with frequency control and reactive power 
control (FCVC). 

The analysis of the impact of wind power on the electromechanical oscillations has been 
done by observing the location of eigenvalues for different wind power penetrations in a 
test system that allows to isolate the oscillations modes. It was showed that as the wind 
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Fig. 10. System response at bus 1 with a wind farm penetration of 20%: (a) frequency, (b) wind farm voltage and 
(c) apparent power of the wind farm. 


power penetration was increased wind farms with FC and with FCVC control strategies 
highly increased the eigenvalues damping for both kinds of oscillations, inter- and intra- 
area ones. Then, these strategies contribute to the damping of the electromechanical 
oscillations. This was corroborated with simulation results. Also, the reactive power 
control allowed to regulate the voltage with a slight decrease in the performance of the 
transient response. Also, simulation results showed that, with a reasonable control effort 
(incremental wind farm apparent power), is possible to highly damp electromechanical 
oscillations. 

Finally, note that, due to the basic characteristic of the control strategies considered, 
it can be expected that investigations on other network configurations will led to 
similar results, counteracting adverse effects of high wind energy penetration in power 
systems. 


Acknowledgment 


This work was supported by UNLP, UNPSJB, CICpBA, CONICET and ANPCyT. 


R.D. Fernández et al. / Renewable and Sustainable Energy Reviews 12 (2008) 1692-1711 1709 
Appendix A. Wind turbine 


The torque and the mechanical power developed by a wind turbine are given by [23] 


2 
Tor 

T, = FP A/A, (16) 
2 

P= Co), (17) 


where p is the air density, r is the radius of the turbine, v is the wind speed, C,(A) is the 
power coefficient, 1 = Q,r/v is the tip speed ratio and Q is the turbine speed. 

Because of Cp is maximum at the optimum tip speed ratio 2= 2o, to extract the 
maximum available power at any wind velocity, the control system should adjust the 
turbine speed to force: 


Agv 
Qi = Qo = ——. (18) 


At this rotational speed the maximum turbine power Pio and the corresponding torque Tio 
result: 


npr (Qr\3 

Pp = (=) Coos (19) 
npr? (r 3 2 

Tio = “2 1. Cpo% > (20) 


being Cpo = Cp(20) the maximum power coefficient. 


Appendix B. Generating units (DFIG) 


Variable speed wind farms powered by DFIG (wound rotor) are the only ones 
considered in this work. The rotor winding is fed using a back-to-back voltage source 
converter as depicted in Fig. 11. 

The operation of an induction machine can be analyzed using the classical theory of 
rotating fields and the well-known d-q model [24]. In this way, the general reference frame 


Transformer 


Bidirectional 
converter 


Fig. 11. Generating unit: double fed induction generator (DFIG) with rotor converter for variable speed 
operation. 
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model for the induction machine is expressed through the following matrix equation: 


Usx R; + pL, —WgL, pL —Wglhm icy 
Usy MgLs Rs + pL, OgLm PL» lex 
Ury = PLi —(Og — o;)Lm R, + pL, —(@, — Or)Lr iy | 
Ury (Mz — Or) Lm PL (Mg — Oy) Ly R, + pL, Irx 


(21) 


where L, and L, are the stator and rotor inductances, Lm is the magnetizing inductance, R, 
and R, are the stator and rotor side resistances, œ, is the electrical rotor speed, œg is the 
reference frame speed, usx and usy are stator voltages and ux and ur are rotor voltages, 
finally isx, isy, ix and iy, are the stator and rotor currents. If œg is equal to the line 
frequency œe, the synchronously rotating reference frame model is obtained. 

The electromagnetic torque Te, the stator active (P.) and reactive (Q.) powers of a 
DFIG are given by 


Te = 3 PL (isxiry = isylrx), me 
Pe = 3 (Usxisx € Usyisy), 
Q: = 3 (Usyisx a Usxisy), 


with P the number of poles. Because of vector control, it is possible to control active and 
reactive powers independently [24,25]. 


Appendix C 


Wind turbine: P = 700kW at v = 11.75 m/s; r = 22m; p = 1.224; Cp(4) = 1.206 x 107547 
—2.71 x 107446 + 2.0399 x 107345 — 6.519 x 107344 + 1.257 x 107723 — 1.54 x 1077274 
2.85 x 10772 — 4.05 x 107. 

DFIG: V = 690 V; R, = 0.0067 Q; P = 4; L, = 0.0075 Hy; Lm = 0.0061 Hy; L, = 0.0084 
Hy; R, = 0.187 Hy; a = 0.3806 turns ratio; Ji4mi = 536096.59 Kgm*; N = 51. 
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